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The viability of achieving gravitational consistent braneworld models in the framework of 

a f(R) theory of gravity is investigated. After a careful generalization of the usual junction 

conditions encompassing the embedding of the 3-brane into a f(R) bulk, we provide a pre- 

Ih ■ scription giving the necessary constraints in order to implement the projected second order 

bJQ' 

effective field equations on the brane. 
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In the last decade braneworld models performed a genuine branch of research in high energy 
physics. In particular, the two papers of Randall and Sundrum envisage our universe as a 

warped brane embedded into a five dimensional spacetime, the bulk. Several extensions and/or 
modifications of this original set up has been proposed in the literature. These extensions include, 
for instance, other dimensions on the brane and in the bulk, the possibility of additional fields in 



Among the possible modifications 



the bulk, and the consideration of modified bulk gravity 
in the bulk gravity, the consideration of higher derivative models was extensively analyzed (for 
comprehensive reviews in f(R) gravity see 6|). In Ref. it is shown that f(R) theories may help 
in avoiding the fine-tuning problem of the Randall- Sundrum [jj model. Besides, several interesting 



cosmological aspects of f(R) braneworld models were investigated It was also demonstrated 



'Electronic address: caramcs@if.uff.br 

^Electronic address: emilia@if.uff.br 

"''Electronic address: hoff@feg.unesp.br;hoff@ift. unesp.br 



2 



recently that in the five-dimensional f{R) braneworld models scope it is possible to obtain a 
model which solves the hiera rchy problem without the necessity of a negative tension brane in the 
compactification scheme 
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When dealing with gravitational aspects of braneworld models, one naturally faces the possibil- 
ity of generalizing the usual field equations on the brane. In fact, even for a bulk respecting solely 
the General Relativity theory, the projected equations on the brane are modified leading to subtle 
but important departures from standard scenarios. 

Our aim in this work is to obtain the appropriate junction conditions and the effective second 
field equations on a 3-brane embedded in a five-dimensional f(R) bulk. At this point, some remarks 
on the previous literature are in the order. The junction conditions for higher-order braneworld 
models with a Gibbons-Hawking term were founded in the Ref. [ijj and the accordingly cos- 
mological equations for the fourth-order f(R) braneworld were obtained in fl3 |. In these precise 
and encompassing papers, as well as in 7|, the junction conditions are obtained by means of the 
conformal transformation relating f(R) and scalar-tensor theories fis| . In fact, there are many 
advantages in the conformal theory since the resulting field equations are again of second order 
(with matter). As it is well known, the conformal transformation performs an isomorphism be- 



tween the space of solutions of the two conformally related theories |16l | , providing a mathematical 
equivalence between them. It is important to remark however that a mathematical equivalence is 
not necessarily a physical equivalence. In fact, if there are matter fields (different from the scalar 
field of the conformal transformation), the conservation equations are different in the two related 



13]. Therefore, 



frames and the advantages of the conformal transformation in this case is doubtful 
as we shall obtain the effective projected equations for the most general case (with matter on the 
brane), we prefer to reobtain the junction conditions via the usual brute force method. Only in 
Sec. Ill we particularize, for simplicity, our equations to the brane vacuum case and therefore we 
shall approach the effective field equations in a conservative and gene ral way. Finally, we point 
out that this program was carried out in Ref. [3]. Nevertheless in 



18| . and subsequently in 



the same junction conditions as in the General Relativity case were used, which is, in fact, only a 
(non rigorous) over simplification, since the junction conditions depends on the embedding gravity 
space. As we are interested in the effective second order (Einstein like) brane gravitational equa- 
tions our procedure towards the suitable junction conditions is well defined. We remark, however, 
that for higher order equations on the brane it is necessary to implement a more powerful method 
to improve the junction conditions (2^. Recently, the necessity of the new junction conditions was 



circumvented in some cosmological applications to the braneworld case 



2I| 



3 



This paper is structured as follows: in the next Section we briefly review the scalaron in the 
bulk. In the section III we obtain the junction conditions for the f(R) bulk gravity without using 
the conformal transformation method. It is shown that in the appropriate limit we recover the 
usual Israel-Darmois equations. Moving forward we obtain the effective second order projected 
field equations on the brane. In Sec. IV we implement several necessary constraints in the bulk 
f(R) gravity and study the correspondent implications on the brane equations in order to set a 
physically acceptable and viable model. Finally, in the conclusion we summarize our main results 
and give a simple example of how to apply the general projection to a specific (toy-)model. 

II. FIELD EQUATIONS AND THE EXTRA SCALAR DEGREE OF FREEDOM 

The gravitational field equations in the metric formalism for a five-dimensional bulk in a f(R) 
gravity is given by 

f'(R)R AB - \gABf{R) + 9A B af(R) - VaVb/'CR) = klT AB , (1) 

where A = 0, ...,4, f'(R) = df(R)/dR and fc| is the five-dimensional gravitational coupling con- 
stant. An important feature arising in f(R) gravity can be observed when we take the trace of the 
above equation, which reads: 

f'(R)R - ^f(R) + 4D/(JR) = k\T . (2) 

We verify that the trace of the field equations has a fully different meaning of that one it has within 
General Relativity theory. In the latter case, the trace of the field equations is a mere algebraic 
relation involving the Ricci scalar R and the trace of the energy-momentum tensor T. On the other 
hand, in the context of a f(R) gravity in the metric formalism, the trace of the corresponding field 
equations, given by the equation @ is a dynamical equation, scalar field-like, with /' = f'(R) 
playing the important role of the scalar field emerging as a consequence of the modified gravity. 
This scalar field is commonly called scalaron. In other words, this is a further scalar degree of 
freedom (propagating across the entire bulk), which the modification of the General Relativity 
theory gives rise. The equation ([2]) can be written as 

□/' - (y - ir) = klT ■ (3) 

Let us notice that such an equation has the usual form of a Klein-Gordon-type equation for a scalar 
field 4>, which is 

u <t>-^r = S, (4) 
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where V corresponds to a potential associated with the scalar field and S is a source term. Let us 
recall that the mass of the scalar field is obtained through the second derivative of the potential 
evaluated at a minimum of the field 

d 2 v{4>) 



(5) 



For the equation (|3|) we have an effective potential, V e tt, associated with the scalaron, defined 



as: 



dV f 



eff 



df 



5/ f'R 



(6) 



As previously performed by 



221 ]. by means of such expression, we can obtain the mass of the 
scalaron, by considering that its minimum lies at high curvature regimes, where it is assumed that 
\f"R\ <C 1. Since General Relativity succeeds to describe correctly the phenomena observed in the 
universe at high redshifts, when the high curvature approximation is important, it is reasonable to 
choose /' = 1 as the minimal value of the scalaron. So, by adopting such approximation and using 



the equations © and © the scalaron mass will be 



m), 



3 J_ 

8f ' 



(7) 



which imposes positiveness for the second derivative of f(R), i.e /" > 0, in order to guarantee a 
tachyons-free theory. This is the first constraint to be satisfied by a viable f(R) in our model. 

III. THE APPROPRIATE JUNCTION CONDITIONS AND THE PROJECTION 

SCHEME 



The junction conditions 'measure' (via its extrinsic curvature) how a co-dimension one hyper- 
surface is embedded into a given bulk. When the embedding space gravity is given by the General 
Relativity theory, the junction conditions are the usual Israel-Darmois matching 



23, 



241 ] . Here we 



shall derive the appropriate conditions for a bulk in the framework of a f(R) gravity. 
The equation ([1]) may be recast into a more familiar form, reading 

Rab - -^9AbR = T h AB \ 



(8) 



where 



nnbulk 
l AB 



klT AB 



^Rf'(R) - \f{R) + Df'(R) I g AB + V A V B f'(R) 



(9) 
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Following the standard approach, we start from the Gauss equation relating the brane curvature 
tensor Rp^§ to its higher-dimensional counterpart 



RBC D h a A Kh c 'hf + KZKps - KfK M , 



(10) 



where is the brane metric, related to the bulk metric by gAB = h>AB + nA«B being tia an 
ortonormal vector in the extra dimension direction and = h^h^hV ' a^b is the extrinsic cur- 
vature of the 3-brane. With the aid of Eqs. (|10p and ([9]) it is simple to obtain the Einstein tensor 
on the brane 



G 



1 AB n n n v + | 1 AB n n 



bulk^AB 



%ulk 



\h 



[IV 



1 



+ KK^ - KZK av -lhp,(K*- K a ?K a 



■E, 



(11) 



where = Cqq^uati hf[h„, the usual projection of the 5D Weyl tensor (Cg CD ) encoding 
genuine extra-dimensional gravitational effects. 

It is conceivable, for the projection scheme purposes, to decompose the bulk stress tensor into 
Tab = —A-9AB + SABOiS(y), where A is the five dimensional cosmological constant, Sab the stress 
tensor on the brane, the Dirac distribution localizes the brane along the extra dimension 1 and a 
is a constant parameter with dimension of length, introduced here in order to allow Sab to have 
the same dimension of Tab- Moreover, the brane energy- momentum tensor is also decomposed in 



a similar fashion, S^ u 



-Xhftv + Tavi m which A is the brane tension and t^ v stands for matter 



on the brane. From Eq. ([9|) and the above decomposition of the stress-tensors it is easy to see that 
the first three terms appearing in the right hand side of Eq. (jlip are 



C — qpbulk^A^B , rffoulk^AB -l 
t,tiv = 1 AB tl a' L u+ 1 AB n n h 



\rjnbulk 



fJ,V 



f'{R) 



2k 2 5 A + (Rf'(R) - f(R)) + 2Df'(R) 



5k 2 5 Af(R) 



in such way that 



V A V B f(R)n A n B 



f'(R) 



V A V B f'(R)rfhZ, 



(12) 



G 



K°K av 



E 



(13) 



So far, the obtained equations agree with Sec. II of Ref. [13] . However, from now on, we shall 
generalize the usual junction conditions in order to accomplish the f{R) bulk gravity and determine 
the extrinsic curvature terms in the projection procedure. In general grounds, we are assuming that 



1 We note, by passing, that the functional form of Tab could also be complemented via the addition of a bulk matter. 
In particular, the junction conditions to be derived in what follows allow one to do so. 
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the normal derivative of Kn U may become large when compared with its variations along the brane 



dimensions 



25f | . Therefore, the relevant discontinuity is computed from [K^ u ] = — K , where 



Kj^ u means the limit of Kn U approaching the brane from the ± side. These considerations can 



be made more precise by the following construction [26|]: consider the brane as a timelike surface 
intersected orthogonally by geodesies. In a given coordinate system it is possible to set the brane at 
y = in agreement with the previous decomposition of the bulk stress tensor. Hence it is possible 
to write ua = QaD- Moreover, it is quite convenient to define another coordinate system, say z^, 
installed upon the brane such that the hypersurface may be parametrized as x A = x A (z^). In this 
vein, it is possible to define the soldering jacobian e A = dx A jdz^ ', tangent to the curves belonging 
to the brane. Note that, by construction, e A iiA = and thus, e A shall act as a projector onto the 
brane directions. Now, it is well known that, with the help of the Heaviside distribution Q(y), the 
bulk metric may be decomposed into different metrics on both sides of the infinitely thin brane 
by 27] gAB = @(?/)5ab + ®(~v)9ab- ^ ne terms g\ B mean the metric in the ± side and the G(y) 
distribution obeys the rules 2 (y) = @(y), @(y)Q(—y) = and d@(y)/dy = 5(y), which makes the 
calculations quite manageable. 

By decomposing the metric via the Heaviside distribution all the relevant geometrical quanti- 
ties may be obtained. For instance, the first derivative of the metric reads qab,c = ®(y)9~AB c 
®(~u)9ab c + < Kz/)[3AB] n Cj leading to the imposition of the Darmois junction condition [gAB] = 0, 
since the product ®(y)5(y) (which would appear in the connection) is not defined in the distribu- 
tional calculus. It is important to note that the Darmois condition guarantees the continuity of its 
tangential derivatives. Therefore the only possible discontinuity of the metric derivative shall be 
along the extra dimension: [gAB, c] = K ABnc, where the kab tensor will not appear in the final 
condition. 

Following this line it is possible to connect the delta part of the Einstein tensor in the brane 



with the delta part of the decomposition f% k = ®{y)f^ B + G(-y)Tj B + 8{y)f A B |27|]. Now, by 
construction, we have T^^e^e^ = T^ u in such a way that the Eq. ([9]) reads 



T 



k 2 5 S^ u - -Rf'(R) - -f(R) + Df'(R) V + e A e*V A VBf(R) 



(14) 



f'(R) 

Note that in this last equation the first term is given by S^ v (the stress tensor on the brane) and 
not by T^ u , since one must respect the constraint T^ u n u = 0. Now, since the geometrical part of 
the Israel matching condition is not modified, the delta part of the Einstein tensor will be the same 



of the standard case [27j. Therefore, the appropriate junction condition reads straightforwardly 

1 



[w\ 



(15) 



with T^y given by Eq. (|14p . Hence, it is possible to express the jump of the extrinsic curvature in 
terms of the brane stress tensor plus corrections coming from the f(R) bulk. It is given by 



a 



[IV J 



f'(R) 

whose trace reads 



1 



1 



ms^ - + aKARfiR) - f(R)] + 4e^ A v B f'(R) 



(16) 



[K] 



a 



k 2 5 S - 2Rf'(R) + 2f(R) - 3D/'(i2) 



(17) 



3f'(R) 

Note that by taking f(R) = R (f f (R) = 1) we recover the usual General Relativity case, as 
expected. 

From Eqs. (|15p and (|16p it is possible to complete the projection procedure, writing down 
Eq. (|13p on the brane. In order to implement it we notice that the Z<i symmetry imposes the 
condition n\ — > —n~l, leading to Kt„ — > —K 



[XV 



[xv 



28] . Therefore, it is possible to calculate the 



extrinsic curvature terms of Eq. (|13p on the brane (henceforth we suppress the label ±, since the 
dependence of in Eq. (|13p is quadratic). The calculation of the extrinsic curvature terms is 
easy but tedious. We will suppress some details and furnish the main results in terms of the brane 
stress tensor. They are 

„2 f l4 r.t.2 



KK 



or 



H v 



12[f'(R)f 



klSS^ - -^-h^S 2 + -^h^S[Rf'{R) - f(R)] + k 2 Se^V A V B f'{R) 



1 



- 2k 2 S^[Rf{R) - f(R)] - -h^[Rf'(R) - f(R)} 2 - 2[Rf'(R) - f(R)]e^V A V B f'(R) 



3k 2 S^nf(R) + k 2 h^SUf{R) - -h^[Rf'(R) - f(R)]nf(R) 



3Uf{R)e^V A V B f{R)\, 



(18) 



k;k uct 



= jfcf - ^SS^ + lh, u S 2 J + M V [i?/'(i?) - f(R)] - ^Sh^[Rf(R) ~ f(R)\ 

+ k 2 (S^ + SZtf)e*V A V B f'(R) - h 2 Se^V A V B f'(R) + 1 h» v [Rf{R) - f(R)} 2 



+ -[Rf'(R) - f(R)]e^V A V B f(R) + e^V c V A f'(R)V c V B f'(R) , 



(19) 
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a 



klS 2 - Ak 2 S[Rf'(R) - f(R)) - 6k 2 SDf'(R) + 12D f'(R)[Rf'(R) - f(R)} 



36/'(i?) ^ 
+ 4[Rf'(R)-f(R)] 2 + 9(nf(R)) 2 



(20) 



and 

K a/3 K, 



Q- 



aB 



\S 2 ) - ^S[Rf'(R) - f(R)} + 2k 2 S^e A e B V A V B f'(R) 



2k 



1 



1 



=fSUf'(R) + -lRf'(R) - f(R)] 2 + ^[Rf'(R) - f(R)]nf(R) 



+ V A V B f'(R)V A V B f'(R) y 



(21) 



Now, from Eqs. (|20p and (|2ip . it is possible to calculate the scalar part of the extrinsic curvature 
terms in Eq. (113|> . Taking into account the decomposition Suu = —Xh^ + t^ u we have 



a 



K2 - K ^ = mR) { 



t 2 2\r AX 2 a 



i 



A' IT 



(r-4X)[Rf'(R)-f(R)] + (Df(R)f 



+ - f(R)) + -[i?/'(i?) - /(i?)] 2 - 2k 2 T^e A e B V A V B f(R) 

- V A V B f'(R)V A V B f'(R)\. 



(22) 



In a similar fashion, the remaining extrinsic curvature terms of Eq. (|13|) are given by 

„2 r / i \ i 



a 



4/'(/2) 



^5 T^Tjju + -TT, 



//// 



fc 5 V 2fc 5 2 A + [#/'(#) - /(i?)] + 3D/'(i?) 



- \e A e B V A V B f\R) (k 2 r + 2fc 2 A + - /(i?)] + 3D/'(i2)) + ^ 

+ ^r[Rf'(R) - f(R)] - k 2 X[Rf'(R) - f(R)) - \[Rf(R) - f(R)] 2 

- k 2 XUf'{R) + k 2 rUf\R) - \{Rf{R) - f{R)}Uf'{R) 

+ kl(r°e A + r^e A )e B V A V B f'(R) + e A e B V c V A f'(R)V c V B f'(R)^. 

The last two equations allow one to bound all the extrinsic curvature terms in Eq. (|13p by 
means of the brane stress tensor components, as well as f{R) and its derivative terms. We shall 
write it explicitly, since these terms encode the junction condition departure from the usual case. 
Hence, organizing similar terms for convenience we have 



(23) 
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K K — K a K 



K 2 - K a/3 K, 



aB 



4/'(i?) I 5 ! M 3 6 2 



hfj, v T a °T a g 



hl^ u (2k$Tp,\ + [Rf'(R) - f(R)] + 3D/'(i2)) - ±e A e B V A V B f'(R) (k 2 r + 2k 2 X 



+ [Rf'(R)-f(R)]+3Df'{R)\ + k 2 (r^ + r^ )e B V A V B f'(R) + e A e B V 



x V A f{R)VcV B f'(R) + 3 V 



2fc 5 4 Ar + K 5 4 A 2 + \k 2 r[Rf{R) - f(R)} - \k 2 \[Rf{R) - f(R)} 



+ \[Rf\R) - f{R)f + [Rf'(R) - f(R)]af(R) - k 2 \Of'(R) + k 2 rUf\R) + \{Uf'{R)) 2 
- 3k 2 r^e A e B V A \/ B f'(R) - *-V A V B f'{R)V A V B f\R) 

Now it is possible to write down the effective gravitational field equation on the brane. Thus 
from eqs. (JT2J) and ([MI), eq. reads 

1 



(24) 



G 



'■^■effhfiu + k^TT^ u + k 2 G e ffT hLV + 



V A V B f(R)h A h B + G eff e A e B V A V B f'(R) 
4[m p (< e ^ + <e^)e^V A V B /'( J R) - _^^V c V A / , (U)VcrVB/ / (JS) 



/V 2 



fii/i 



(25) 



where 



and 
A 



G 



a 



2k 2 X + [Rf'(R) - f(R)] + 3Df'(R) 



TT 



( r CT T + -TT + T 2 --/i T afi T « I 

4[/'(i?)] 2 I M 3 6 2 ^ / 



(26) 



(27) 



e// = 4[//(jR)]2 S pfAr + ^A 2 + -k 2 r[Rf(R) - f(R)] - -k 2 X[Rf'(R) - f(R)] + -[Rf'(R) - f(R)] 2 
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1, 



- M A D/'(fl) + ^rD/'(i?) + ^(□/( J R)) 2 - ftf-r^ef V, 



Z.2 

Kg 

3 



-V A V B f(R)V A V B f'(R) } + 2fc 2 A + [i2/'(2J) - /(#)] + 2D/'(i2) 



5fcfA/'(i2) 



VAV B /'(i?)/n B 



(28) 



In the next Section we shall implement a few important conditions coming from braneworld 
models. These constraints acts supplementing the covariant approach of this Section, leading to 
well defined five-dimensional braneworld models within f(R) gravity. By now we just remark, by 
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passing, an important output encoded in the Eq. (|26p . Note that the f(R) bulk gravity is felt on 
the Newtonian effective constant. Thus, in order to guarantee a positive G e // it is necessary that 2 

2k 2 5 X + [Rf(R) - f(R)] + 3Df(R) > 0. (29) 

Therefore Eq. (|30p is the first nontrivial constraint which must be respected by a gravitationally 
viable braneworld model in the f{R) bulk context. 

IV. CONSTRAINING AND ENABLING f(R) PROJECTED BRANEWORLDS 

In order to better appreciate the relevant constraints which have to be respected for a theoret- 
ically and physically interesting model, we shall look at the simplest case of the vacuum on the 
brane, e. g. t^ v = 0. In this case, equations ([25]) and ([28]) give 

= -k e}f h^ + j^V A W B f\R)hffi + G eff e^V A V B f'{R) 

^e^ c V A f'(R)V c ^Bf'(.R) ~ E^, (30) 



with 



*[/'(*)] 



° 2 I 1'4\2 1 ;„2 \\ r> ft I n\ «n\l i 1 [nf/D\ ;/n\l2 ^ 



A e// = vjf^I | 3 fe f A " g*5 W(*) - /(*)] + ^[^/'(i?) - W + ^[Rf'(R) ~ f(R)]Of'(R) 
- fxaf(R) + l -{Uf'{R)f - l -V A V B f (R)V A V B f(R)\ + 2k%A + [Rf'(R) - f(R)] + 2D/'(i?) 



5fc 5 2 A/'(E), 



VAVs/'CRKn* (31) 

In general, when dealing with /(i?) gravity based braneworld models, there are two types of 
physical inputs which may be accomplished by a well defined model: on the one hand, there are 
constraints coming from braneworld theory (such as the mentioned positivity of the brane tension) 
and, on the other hand, there are additional conditions arising from the requirement of a viable f(R) 
gravity model, necessary to fit experimental and observational data. Obviously, the implementation 
of the former type of conditions is indeed necessary. However, the direct use of the f(R) constraints 
in our case is a naive approach, since these constraints are valid for usual four-dimensional f(R) 
cosmological models and we are dealing with f(R) gravity in the bulk. Besides, the obtained field 
equation ( (|3"0"]h for instance) is quite different from the standard case. In this vein, our approach 



2 The brane tension A is the proper vacuum energy for isotropic branes. This reinforces the fact that A is a positive 
parameter 



11 



here is to follow the clue of the braneworld constraints and complement it with some additional 
conditions appearing in the f{R) theory, without making reference to any particular model. 

In the final part of the Sec. II we mentioned that the brane tension should be positive; also, in 
order to get a attractive gravity on the brane we must impose that G e ff > 0. We can use equation 
([2]) to remove df'(R) and express (f30|) solely in terms of some constants, f(R) and its derivatives 
with respect to R. In doing so we have 

^/' + £/ + ^ 5 2 |A|+2fc|A>0. (32) 

In a pursuit for a physically interesting braneworld model, it is conceivable to insist in having 
an AdS bulk. Therefore one shall write A — > — |A|. Besides, the minuteness of A e ff (it shall be 
lower than 10 -120 ) is well known and may be associated to an important constraint. 

Let us follow a similar procedure in (13ip by using (JTJ) and (JSj) in order to remove the terms of 
derivative of f(R) with respect of the coordinates, i.e \V A V B f'(R)V A V B f'(R), V A V B f'(R)n A n B 
and n\f'(R) ending up with an equation depending only on f(R) and its derivatives. With such an 
upper bound in the Eq. (|3ip the resulting constraint reads 

\kt\ - ±kl\Rf + ^kt\f + I/ 2 + hi\A\fR + A/jg|A| + lf 2 R 2 + |*g|A| 2 + ^/% 2 |A| 

+ ~ 2 f' 3 R + l^f- ^!A|A| - \f' 2 RABR AB + - ^-R AB n A n B = . (33) 

or 2 or 12 2 or or 

The point to be stressed here is that any five-dimensional braneworld model based upon an AdS 
f(R) bulk must respect Eqs. (I32p and (J33J) in order to guarantee a well defined model from the 
gravitational point of view. Of course, Eqs. f|32|) and (|33j) are rather nontrivial and it is not expected 
that they could be satisfied by chance. We remark, by passing, that in the General Relativity limit 
none of the constraints imposes any bound on the curvature scalar, as expected. 

V. FINAL REMARKS 

In this paper we explored the viability of achieving a gravitational consistent braneworld scenario 
in the framework of a f(R) theory. For this purpose, we worked in the most general context without 
considering any particular model in order to derive the appropriate junction conditions to obtain 
a prescription on how to project the effective second order equations on a 3-brane embedded in a 
five-dimensional bulk. Besides providing a rigorous derivation of the junction conditions, our main 
result here is that any five-dimensional braneworld model based upon an AdS f(R) bulk must 
respect Eqs. (|32p and (|33p in order to guarantee a well defined and gravitationally viable model. 
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Finally, let us give a brief exampliflcation of how to apply equations (|32p and (|33|) to a particular 
(toy-)model. Let us consider a case in which we are dealing with a theory of gravity consisting of 
a small modification on General Relativity: f(R) = R + <p(R), <p(R) *C R It is reasonable to single 
out models obeying this condition, since it is expected for an arbitrary f(R) not to be so much 



2J. 



different of General Relativity, which ensures the fulfilment of the viability conditions listed in 
By adopting such an assumption we will be led to the approximations below 

f = R2 ( 1 + 2 i)i (34) 

f' 2 = 1 + 2^, (35) 

f' 3 = 1 + 3<p' (36) 

and 

f 2 f = i?(l + 2^ + |). (37) 



So, we will get a linear form for (|33p which reads 



1 R AB R AB ~ -^R AB n A n B = . (38) 



ai R 2 + a 2 R + a 3 - R AB R AB - ^R A Bn A n B + (6i + 6 2 i?) ^ + ci^ 2 + c 2j R + c 3 

cr / a/t 



2 a 
With the constants ai, a 2 , 03, 61, 6 2 , c\, c 2 and C3 given by 

ai = \, a2 = \ k ^~ je k ^ X + ^ ' a3 = 5 |A| ' 61 "ll' 62 S 72^ A+ 24^' A ' + 

ci = - , c 9 = -£|A - -Ag|A| + , c 3 = -fc 2 A + -fcf |A| 2 + Sa| - -fc§A|A| . 
32 ' 72 5 8 51 1 2a 2 3 5 8 51 1 a 2 1 1 12 5 1 1 

The equation (|38[) can just be solved in the case in which the scalars RabR AB and RAB"n A n B are 
expressed as functions of the Ricci scalar R. Such equation represents an important condition of 
viability for models possessing the functional form f(R) = R + <p(R) and obeying the constraint 
<p(R) <C R. Apart of that, it is quite conceivable to implement a typical warped line element, ex- 
pressing RabR AB an d RAB nAnB i n terms of the warp factor. In this vein, the resulting differential 
equations for the warp factor in the coefficient of dipjdR and in the second term of Eq. (|38|) may 
be used to classify warped spaces accordingly to the correction ip. 
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